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. We apply power series expansion to symmetric multi-well oscillators bounded 

by two infinite walls. The spectrum and expectation values obtained are compared 
with available exact and approximate values for the unbounded ones. It is shown 
that the method is capable of producing to a high accuracy the eigenvalues, eigen- 
functions, and the expectation values (x 2k ) of the corresponding unbounded ones 
£h ' as the separation between the two infinite walls becomes large. 
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1 Introduction 



Potential energy functions with degenerate classical minima separated by potential bar- 
riers are frequently encountered in many areas of physics, such as molecular physics, 
nonlinear physics, field theory, and cosmology. A widely investigated model, by variant 
methods ([1]— [19]), is a one-dimensional symmetric double- well potential, in which the 
degeneracy of the low-lying doublets is broken by quantum tunneling effects [20]. 

Recently, calculation of energy splitting of a one-dimensional symmetric multi-well 
potential energy functions has received a considerable interest ([21]-[27]). This interest 
stems from the fact that the topology of the multi-well potential is completely different 
from the double-well case. Each well, except the two outer ones, is surrounded by two 
barriers. In the case of odd number of wells, the positions of the minima of the potential 
are connected by the spatial inversion symmetry of the potential that leaves the position 
of the minimum at the origin invariant, for even number of wells ,however, these minima 
are connected by the inversion symmetry. These topological properties are relevant for 
calculations based on path integral methods using instanton solutions[28, 29, 30]. 

In an earlier paper [19], we have applied the traditional power series expansion method 
with boundary conditions in the finite interval to one- dimensional anharmonic potential 
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functions having single and double-well. We have shown that the method is capable, even 
for a reasonable size of the interval, of producing to a high accuracy the eigenvalues and 
eigenfunctions of the corresponding unbounded ones. In particular for the deep double- 
well case, with potential function (V(x) = —fi 2 x 2 + x 4 ), we have shown that the splitting 
of the energy levels of the ground and first excited states becomes apparent only after 26 
significant digits for the mass parameter fj, 2 = 25, and for fi 2 = 35 the splitting shows up 
after 42 significant digits, while for /i 2 = 50 the splitting starts after 72 significant digits. 

In the present paper we extend our previous calculation to multi-well potentials and 
compare with available exact and approximate values recently obtained by other methods. 

The rest of the paper is organized as follows. In sec. 2 we apply the method of power 
series expansion to anharmonic multi-well potentials, and present our calculation of the 
eigenvalues, eigenfunctions, and expectation values of even powers of the x-coordinate and 
compare with that recently obtained by the improved Hill determinant method [21, 23], 
and with the exact values calculated by supersymmetric quantum mechanics [21]. It is 
shown that the method yields extremely very close results to the exact ones. Finally in 
sec. 3 we give our conclusion. 



2 Calculation of energies and expectation values 

In order to calculate the eigenvalues and eigenfunctions of the anharmonic multi-well 
potential bounded by infinitely high potential walls located at x — ±L, one needs to solve 
the eigenvalue equation (in units h = 1 , 2m = 1) : 



d 



'2 



dx< 



+ E- V(x) 



*(x) = 0, (1) 



with the boundary conditions ^(±L) = 0. In the present work we consider potentials 
V(x) of the form [22, 23, 27]: 

TV 

V(x) = J2b2kX 2k , (2) 
fc=i 

Here the coefficients (62m) are real with & 2 tv being positive, and N = 3,4, 5. 
Making use of the power series expansion 

00 

= ]T a n x n (3) 

n=0 

in eq. (1) we obtain the following recurrence formula for the expansion coefficients: 
— E a n _ 2 + i>2 a n _ 4 + 64 a n _ 6 + 6 6 a n _ 8 + b$ a n _i + 610 a n -i2 , „ -. 

On= 7 7T , n^0,l 

n [n — 1) 



a n = 0, n < 



(4) 



The symmetry of eq. (1) divides the solutions into two types even and odd. The even 
solutions can be obtained by imposing (ignoring normalization) ao = l,a\ = 0, while 
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the odd ones by imposing ao = 0, ai = 1. The energy eigenvalues (i?) are obtained from 
the condition ^(L) = in both cases. Since we are dealing with potentials admitting 
power series expansion for \x\ < L, the power series solutions of *&(x) are, according 
to a well known theorem in differential equations, convergent [31]. Here we present our 
calculation of the energy levels, wave functions and expectation values for the potentials 
given in eq. (2), using eqs.(3),(4) and the boundary condition at x = L. 

Table 1 shows the first four energy eigenvalues of three-, four-, and five- well oscillators 
for the bounded potential-wells as compared with the obtained values reported in [22, 23] 
for the unbounded case. The obtained values are extremely very close to the exact ones. 

Table 2 contains a comparison of the expectation values (x 2m ),(m = 1---5) with 
the reported ones and the exact values obtained from the corresponding supersymmetric 
wave functions. For completeness we also present in Fig. 1 (A, C) the ground state 
wave functions for the potentials A(x) and C(x) respectively, and (B,D) are for the 
first excited state wave functions for the potentials B(x) and D(x) respectively. These 
coincide with the corresponding exact wave functions (ignoring normalization) obtained 
using supersymmetric quantum mechanics [21]: 




Figure 1: These graphs show the ground state wave functions for the potentials A(x) and 
C(x), also showing the first excited state wave functions corresponding to the potentials 
B(x) and D(x) for the bounded case 
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Potential 


No. of 
wells 


L 


Our results 


Results of 
supersymmetric 
quantum mechanics 


A(x) 


3 


4 


-2.000 000 000 000 000 000 
-1.772 726 698 991 350 330 
2.078 279 891 768 595 361 
5.604 028 342 382 013 654 


-2 


B{x) 


3 


4 


-9.001 720 238 527 719 715 
-9.000 000 000 000 000 000 
0.639 394 262 865 333 280 
1.936 629 224 926 380 607 


-9 


C(x) 


3 


3 


0.375 000 000 000 000 000 
2.357 398 881 839 175 696 
6.988 755 014 467 723 300 
13.880 516 236 71388 190 


0.375 


D(x) 


3 


3 


-0.195 122 059 734 627 597 
1.125 000 000 000 000 000 
5.646 143 524 135 629 302 
12.384 746 745 98 872 574 


1.125 


F(x) 


4 


4 


-0.223 991 055 384 171 854 
0.083 481 557 863 966 793 
1.526 487 708 073 844 797 
3.971 174 256 474 939 459 




G(x) 


4 


3 


-0.096 291 946 230 649 098 
0.672 993 242 745 446 704 
3.111 022 328 724 771 653 
7.038 082 659 880 398 654 




H(x) 


5 


4 


0.807 741 647 209 432 443 
3.277 946 311 571 061 982 
7.667 480 496 116 480 534 
13.578 984 131 990 285 801 





Table 1: The first four energy eigenvalues of three-, four-, and five- well oscillators: A{x) = 
x 2 - 4x 4 + x 6 , B(x) = Ax 2 - 6x 4 + x 6 , C(x) = ^x 2 - fx 4 + x 6 - x 8 + x 10 , D(x) = ^§x 2 - 
f x 4 + x 6 -x 8 + x 10 , F(x) = -x 2 + 2x 4 -0.9x 6 + 0.1x 8 , G(x) = -x 2 + 3x 4 -2x 6 + 0.1x 10 , H(x) = 
2x 2 — 7.5x 4 + 5.5x 6 — 0.877x 8 + 0.04x 10 as obtained in the present work for the bounded cases 
(2L is the width of the well). The underlined numbers correspond to the unbounded ones [23]. 
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Potential 


Expectation 




results 




and eigenstate 


value 


Our results 


in [22] 


Exact 


A(x),V (x) 


(x 2 ) 


1.7042723043 


1.704 


1.7043 




\ x ) 


3.9085446087 


3.908 


3.9085 




(x 6 ) 


10.373497673 


10.349 


10.3735 




(x 8 ) 


30.518356869 


30.212 


30.5184 




(x 10 ) 


97.343955598 


93.210 


97.3440 


5(x),*i(x) 


(x 2 ) 


3.1795525642 


3.179 


3.1796 




(x 4 ) 


11.0386576927 


11.035 


11.0387 




(x 6 ) 


41.0648544887 


40.985 


41.0649 




(x 8 ) 


161.8298653908 


159.909 


161.8299 




(x 10 ) 


670.2814413720 


625.164 


670.82814 




(x 2 ) 


0.45832470069 


0.458 


0.4583 




(x 4 ) 


0.43854420934 


0.439 


0.4385 




(x 6 ) 


0.54740034191 


0.547 


0.5474 




(x 8 ) 


0.79673314349 


0.796 


0.7967 




(x 10 ) 


1.28945196690 


1.288 


1.2895 




(x 2 ) 


0.956841751448 


0.957 


0.9568 




(x 4 ) 


1.194350514116 


1.194 


1.1944 




(x 6 ) 


1.738359600265 


1.738 


1.7384 




(x 8 ) 


2.8134027359576611 


2.812 


2.8134 




(x 10 ) 


4.935043317278 


4.930 


4.9350 



Table 2: Comparison of the expectation values (x 2m ) of the supersymmetric potential 
A(x), B(x), C(x), and.D(x) obtained in the present work for the bounded cases and that of 
ref. [22] for the unbounded ones along with the exact values obtained using supersymmetric 
quantum mechanics [21, 22] 
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*o(a;) 



exp (x 2 — x 4 /4) for A(x) 

exp(-3x 2 /16 + x 2 /8-x 6 /6) for C(x) 

x exp (3x 2 /2 - x 4 /4) for B(x) 

x exp (-3x 2 /lQ + x 4 /8- x 6 /Q) forD(x) 



These calculations show that even for a reasonable size of the interval the present method, 
beside being simple and direct is highly effective. 

3 Conclusion 

In this paper we have applied the method of power series expansion to a variant of 
one-dimensional multi-well potential functions bounded by two infinite walls. We have 
compared our calculation of the low-lying energy levels, wave functions, and expectation 
values with that obtained using improved Hill determinant and supersymmetric quantum 
mechanics. For all these potentials we have obtained results that are extremely very 
close to the exact ones. This is because the power series expansion in the finite range 
is convergent for potentials admitting power series expansion. It is easy to extend the 
present method to other multi-well polynomial potential functions with higher degrees 
and more degenerate classical minima. 
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